Ab initio calculation of low-energy collective charge-density excitations in MgB 2 
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We present ab-initio time-dependent density-functional theory calculation results for low-energy 
collective electron excitations in MgB 2 . The existence of a long-lived collective excitation corre- 
sponding to coherent charge density fluctuations between the boron a- and 7T- bands (air mode) is 
demonstrated. This mode has a sine-like oscillating dispersion for energies below 0.5 eV. At even 
lower energy we find another collective mode (aa mode). We show the strong impact of local- field 
effects on dielectric functions in MgB2. These effects account for the long q-range behavior of the 
modes. We discuss the physics that these collective excitations bring to the energy region typical 
for lattice vibrations. 

PACS numbers: 74.25.Kc,71.45.Gm,73.43.Lp,74.70.Ad 



After several years of intense experimental and the- 
oretical study of the origin of unusually high temper- 
ature superconductivity^ in MgB2, nowadays there is 
a general belief that the superconductivity in MgB2 is 
phonon-mediated with multiple gaps and an exceptional 
role played by E2 g phonon mode in strong coupling to 
electrons in the tubular-shaped a bands at the Fermi 
level, Ep. However, despite the tremendous progress 
achieved in understanding of superconductivity in MgB2, 
many of its fundamental properties related to phonons 
and electron-phonon (e — ph) coupling are still puzzling 
and remain under intense debate. Thus such important 
issues as, e.g., the reduced isotope effect )2j2i£ the nature 
of the strong broad continuum in Raman spectra, the 
unusually broad linewidth of the Ei g phonon mode seen 
in Raman spectra^i the inconsistency between the re- 
sults of Raman scattering and x-ray measuremcnts^^^ 
and the role of non-adiabatic effects ^ 11 ' 1 ^ 13 ' 14 have not 
yet found satisfactory explanations. 

As for the electronic part of the e — ph picture of super- 
conductivity in MgB2, it has received substantially less 
attention than the phonon part. In particular, in all eval- 
uations of the e—ph interaction the adiabatic approxima- 
tion has been used, and even in ab initio calculations^ 
of T c the dynamical Coulomb interaction has been con- 
sidered in its static form. Here we present a detailed 
ab initio study of the low-energy dynamical electronic 
properties of MgB2. We demonstrate that the strongly 
anisotropic electronic structure of MgB2 characterized by 
boron quasi two-dimensional a and three-dimensional tt 
band a 16 i leads to remarkable low-energy dielectric re- 
sponse in this compound: collective modes with a pecu- 
liar sine-like oscillating dispersion appear in the 0-0.5 eV 
energy range. This brings interesting physics to the en- 
ergy region which was thought to be entirely dominated 
by lattice vibrations. 

Our approach is based on time-dependent density func- 



tional theory^ where the non-local dynamical density- 
response function, x, determines an induced charge den- 
sity, p, in the electronic system caused by an external 
potential, u cx t, according to 



p(r,t) = / X (r,t;r',t')v e ^(r',t')dv'dt'. 



(1) 



X is obtained from the integral equation x — X° + 
X°{v + K xc )x, where x°( r J t] r ': O is the response func- 
tion of non-interacting electron system, v(r — r') is 
the bare Coulomb potential, and K xc accounts for 
dynamical exchange-correlation effects. In terms of 
Fourier-transformed quantities the imaginary part of 
the matrix Xgc (Ij w ) can be- evaluated using equation 

^GG'(q. w ) = -i s g n (w)Im[xGG'(q J a; )]r awitntlie spec- 
tral function 5 GG /(q,o;) defined as 
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and Re[x GG / (q, cj)] is evaluated via the Kramers-Kronig 
relation using an energy cutoff of 50 eV. Here the G's 
are the reciprocal lattice vectors, n,n' are band indices, 
the wave vectors k and q are in the first Brillouin zone 
(BZ), the factor 2 accounts for the spin, and f2 is the nor- 
malization volume. In practice, we replace the ^-function 



in Eq. ([2|) by a Gaussian 



v^7 



h 



with a very small 



broadening parameter 7 = 10 mcV. The sum over k was 
performed on a 108 x 108 x 360 grid. The use of such 
a fine k-mesh sampling is crucial to achieve convergence 
of dielectric properties in MgB2 at low energies. The 
one-particle energies e ra k and wave functions ipnk & r e ob- 
tained as self-consistent solutions of Kohn-Sham equa- 
tions using norm-conserving pscudopotcntials 2 ^ and the 
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FIG. 1: (color). Calculated — Imfe^ 1 (q c , ui)] versus uj and q c 
in two energy regions. Calculations include local-field effects 
and the RPA kernel. In (a) the circles and squares mark 
the energy-loss peak positions measured in x-ray scattering 
experiments of Refs. [24] and [25], respectively. In (b) the dis- 
persion of the upper sharp an plasmon peak is described by 
^o-ir = Vctv ■ |sin(|g c )| (green dashed line) with Vmv = 0.48 
eV. The lower broad feature with similar sine-like dispersion 
corresponds to the aa plasmon. Yellow solid line presents 
acoustic plasmon dispersion according to Ref. 
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FIG. 2: (color). Real (red lines) and imaginary (blue lines) 
parts of the dielectric function, e(q c ,ui), and the corresponding 
loss function, — Im[e _1 (g c , lo)], (black lines) for (a) q c = 0.055 
a.u. -1 and (b) q c = 0.322 a.u. -1 evaluated with the use of 
the RPA kernel. Thick (thin) lines stand for results with 
(without) inclusion of local-field effects. In (a) the dashed 
lines are the free-electron gas results for r s — 1.8 a.u. In (b) 
thin dotted lines show contributions to Ime due to intraband 
transitions within the o\ (green) and <T2 (purple) bands. 



exchange-correlation potential of Ref. Hi]. In order to 
elucidate the role of exchange-correlation effects in the 
density-response function, we performed calculations of 
X using two forms of the many-body kernel K xc , namely 
the random-phase approximation (RPA) (A' xc =0), and 
an adiabatic extension of the local-density approxima- 
tion (TDLDA).— 

The calculated loss function, — Im[e _1 (q,a;)]^ di- 
rectly probed in inelastic scattering experiments is pre- 
sented in Fig. [T]as a function of the momentum transfer 
q along the c* axis. In the upper energy range (Fig. 
Hk) it is dominated by a well-defined collective mode dis- 
persing in the 2.5-4.5 eV energy range in excellent agree- 
ment with x-ray experiments! 24 ' 25 The existence of this 
mode for momenta in the first BZ was demonstrated 
in Refs. [26] and H3, whereas its cosine-like oscillating 
dispersion in higher BZ's was recently discovered in a 
joint experimental-theoretical studyj^ 5 . This mode pro- 
duces a strong impact on electrodynamical and optical 
propertie a 28 ' 29 ' 30 of MgB2, however, it has no relevance 
to superconductivity as it affects neither the dynamically 
screened Coulomb interaction at low energies nor phonon 
dispersion i21 

FigureQJ) shows the calculated loss function in the low- 
energy domain that is our main finding. At all q c it is 
dominated by a sharp "<77r" peak with a sine-like disper- 
sion. Figure [2] shows the dielectric and loss functions at 
q c = 0.055 a.u." 1 and q c = 0.322 a.u. -1 . One can see 
how the presence of two types of carriers (in the a and 7r 
bands) around Ep (see Fig. [5^) with a large difference in 
the perpendicular component of the group velocity, v 1 , 
(compare maximal v 1 in the it bands (Fig. [3Jd) with that 
in the a bands (Figs. [3J: and(3]i)) produces in Ime a struc- 
ture consisting of two main peaks (Fig. H^). Thus, while 
the faster tt carriers give rise to a broad structure from 



to 0.68 eV in Fig. [2^ with the main peak at the upper- 
energy side, the a carriers which are moving more slowly 
in the c* direction produce a sharp peak in the low-energy 
part of Ime. The reason of this sharpness resides in the 
fact that within both a bands, the number of states with 
maximal v is greatly enhanced as seen in Figs. [3J;-[3Ji. 
Therefore, the number of intra-band transitions involv- 
ing these fast states is large which strongly enhances the 
a peak in Ime at the higher energies. In turn, this causes 
a dramatic drop in Ree at nearly the same energy which, 
combined with the presence of a local minimum around 
u> = 0.1 eV in Ime, produces a well-defined sharp peak 
in — Im[e _1 ] corresponding to charge density fluctuations 
between the a and tt bands, a cm mode. Despite the pres- 
ence of a non- vanishing value of Ime at the u> where the 
<77r peak appears, its intrinsic width is extremely small, 
being in the meV range (which corresponds to a lifetime 
of few hundreds of femtoseconds), and is almost entirely 
determined by an extrinsic broadening parameter 7. 

Since at small q c the two main peaks in Ime disperse 
linearly with momentum according to their maximal per- 
pendicular Fermi velocity components, Ime cx Vp^J^ -q c , 
the dispersion of the air mode (linked to the upper side 
of the continuum for electron-hole excitations within the 
a bands) is also linear in q c and its group velocity v a7T is 
close to v^jj,^ This corresponds to the acoustic plas- 
mon proposed to exist in MgB2 on the basis of a model 
tight-binding calculation^ 2 - The concept of an acoustic 
plasmon goes back to Pines^ who showed that it can 
occur in a two-component electron plasma consisting of 
slow and fast carriers. The fast carriers (within the it 
bands in MgB 2 ) can act to screen the repulsion between 
slow carriers (within the a bands in MgB 2 ) resulting 
in the appearance of a plasmon mode with a peculiar 
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FIG. 3: (color), (a) Energy band structure of MgB2 along 
some symmetry directions of the Brillouin zone in the vicin- 
ity of the Fermi level, Ef- Maps of the density of states in 
the tv (b), o 1 (c), and 02 (d) bands versus their energy and 
group velocity component along the c* axis. Note the different 
velocity range of (b) in comparison with (c) and (d). 



sound-like dispersion. After the recognition of the role 
in the acoustic plasmon as a possible mechanism for su- 
perconductivity in transition metals^ it has been regu- 
larly evoked for explaining superconductivity in materials 
with unusually high T c . Nevertheless, up to now the is- 
sue remains controversial since the very existence of such 
a collective mode in metals has been confirmed neither 
experimentally nor by ab initio calculations. Only re- 
cently an acoustic-like plasmon was observed in the elec- 
tron energy-loss measurements at a metal surface in ex- 
cellent agreement with the ab initio prediction^ thus 
greatly increasing our confidence in the present results. 

With increasing q c the a peak in Ime, and consequently 
the an peak in — Im[e ], starts to split into two peaks 
(Figs. [TJd and [2p) since for holes, v^££ exceeds 
by more than 20% (Figs. [3Jc) and[3]Jd)). In this case 
the higher (lowcr)-energy plasmon peak in — Im[e _1 ] cor- 
responds to the charge fluctuations between <j\ and n 
(a 1 and 02) bands. Whereas the an plasmon continues 
to be long-lived, the lower-energy aa plasmon has a sig- 
nificantly shorter lifetime due to more efficient scattering 
within the o\ band. We expect the two separate modes 
to exist at any q c , however, when the calculated broaden- 
ing 7 exceeds the energy difference between these modes 
only a single an peak arises in — Im[e _1 ] (Fig. 

In Fig. [TJd one can see how with increasing q c both 
of these plasmon modes reach maximum energy at the 
A point and, as momentum increases further, their dis- 
persions change from positive to negative and approach 
lo = at the T point in the second BZ, in striking con- 
trast to results of Ref. [32|. The origin of this periodic 
behavior resides in the fact that strong local-field effects 
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FIG. 4: (color online). The low-energy part of the bare an 
plasmon (red solid line) along with the calculated^ bare -E29- 
phonon (dashed line) and measured peak positions (open^ and 
filled^ circles). The position of the Raman spectrum peak is 
shown by an orange filled circle at q c = 0. Dispersions of the 
hybrid "an — -&2 9 " modes, P® and P e , are shown by green 
and blue thick dotted lines, respectively. In the evaluation of 
these curves we use a ll an — £29" coupling parameter A=17.3 
meV in order to place the position of the P® mode at q c = 
in agreement with the Raman measurements.—^ 



in MgB2 feed the strength from the small q c modes into 
the charge density fluctuations at large (q c + G z ) in a 
similar fashion as occurs in the case of the higher en- 
ergy mode^ These sine-like dispersing modes continue 
to have strength at momenta in subsequent BZ's and this 
is a direct consequence of a layered MgB 2 structure. Ad- 
ditionally, as seen in Fig. [21 the local-field effects lead 
to a blue shift of the aa and an energy and make the 
an peak more narrow as it occurs for ui corresponding 
to smaller Ime. We also analyzed the role of exchange- 
correlation effects beyond the RPA and found that the 
RPA picture is essentially sufficient for the description 
of low-energy collective excitations in MgB2 . The inclu- 
sion of the TDLDA kernel in the calculation of x{Qc>w) 
leads to only a small (few percent) downward shift of the 
plasmon mode dispersions with a slight reduction of their 
lifetimes, i.e., in part compensating the local-field effects. 

The an plasmon can have deep impact on both low- 
energy-electron and phonon dynamics in MgB2- Our cal- 
culation reveals dramatic modification of the dynamical 
Coulomb interaction in the energy range vital for super- 
conductivity. In Fig. [2] one can see that in the neighbor- 
hood of the an plasmon energy the calculated Rc[e((7 c , ui)] 
differs dramatically from static e{q c ,uj = 0) as well as 
from the free electron gas result. In particular, in a cer- 
tain momentum-energy phase space region, Ree is even 
negative leading to overall reduction of Coulomb repul- 
sion between carriers. This could explain in part the 
observed reduced isotope effect in MgB2^^ It would 
be of great interest to quantify this effect within the ab 
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initio approach^ Note, that the aw plasmon does not 
influence inter-band a-ir and intra-band a scattering due 
to the phase-space restrictions, although could have some 
effect in the intra-band n one. 

The (J7r plasmon can dramatically affect the dispersion 
of the optical phonon modes in the "pathological" re- 
gion of small offering an unexpected explanation of 
widely-discussed discrepancies between the Raman and 
x-ray measurements of the phonon structure in MgB2. 
The ~ 77 meV peak observed in Raman-scattering ex- 
periments and commonly attributed to the boron _E 2g 
phonon mode is always strongly renormalized in com- 
parison with a bare phonon dispersion, whereas the x-ray 
experiments performed for large q c do not see any appre- 
ciable rcnormalization. In Fig. 2] we show the calculated 
dispersion of the bare air mode along with the bare -E 2g 
phonon mode^i compared to the experimental data.— ^ 
One can see how the two bare curves cross each other 
close to the T point. Due to its localization in the boron 
plane the air plasmon strongly interacts with the boron 
Ei g phonon mode resulting in strong hybridization of 
these boson modes. The result shown by green and blue 
dashed lines in Fig. ^demonstrates that the mode seen in 
the Raman experiments is indeed a strongly hybridized 
"cr7r-i?2g" P® mode. Note, that the finite electron life- 
time effec t 11 ' 14 can also contribute to the upward shift 
of this mode. Additionally, a steep air plasmon disper- 



sion at small momenta might explain the existence of the 
strong unstructured background (commonly mentioned 
as being of unknown electronic origin) observed in the 
Raman experiments. This is corroborated by the fact 
that, e.g., in AIB2, where we do not expect the existence 
of such a plasmon, this background does not appear 

In conclusion, our detailed ab initio calculation demon- 
strates the existence in MgB2 in the 0-0.5 eV energy 
range of hitherto unknown long-lived collective mode 
that corresponds to coherent charge fluctuations between 
the boron a— and 7r— bands (an mode) with striking pe- 
riodic sine-like dispersion. This mode shows an acoustic- 
like behavior at small momenta where it can strongly 
interact with optical phonons. Additionally we find at 
slightly lower energy a more strongly damped mode that 
corresponds to charge fluctuations between two different 
a bands (a a mode). Both these modes have profound 
impact on the low-energy dynamical Coulomb interac- 
tion which should be explicitly taken into account in ab 
initio theories of superconductivity to have the predictive 
power. 
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